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We report the results of a study of charged B decays to the D^tt^tt^ and D'^tt^tt^ final states 
using complete D'*-* reconstruction. The contributions of two-body B D**tv decays with narrow 
(j=3/2) and broad (j = l/2) D** states have been determined and the masses and widths of four D" 
states have been measured. This is the first observation of the broad Dq" and D'l mesons. The 
analysis is based on a data sample of 65 million BB pairs collected in the Belle experiment. 

PACS numbers: 13.25.Hw, 14.40Lb, 14.40.Nd 



INTRODUCTION 



B decays to Dt: and D*t: final states are two of the dominant hadronic B decay modes and have been measured 
quite well 0. In this paper we study the production of I?- meson excited states, collectively refered to as D**'s, that 
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are P-wave excitations of quark-antiquark systems containing one charmed and one light (m, d) quark. The results 
provide tests of Heavy Quark Effective Theory (HQET) and QCD sum rules. Figure shows the spectroscopy of 
D-meson excitations. In the heavy quark limit, the heavy quark spin Sc decouples from the other degrees of freedom 
and the total angular momentum of the light quark jq = L + Sq is a good quantum number. There are four P-wave 
states with the following spin-parity and light quark angular momenta: 0+(jq = 1/2), = 1/2), l^(jg — 3/2) 

and 2+(jq = 3/2), which are usually labeled as Dg, DJ., Di and D2, respectively. The two jq = 3/2 states are narrow 




FIG. 1: Spectroscopy of D-meson excitations. The lines show possible single pion transitions. 



with widths of order 20 MeV and have alreadybeen observed [3 S IJ, ISl la l2l la ISl, llOl llH UM ■ The measured values 
of their masses agree with model predictions |l3l ll-d llSl llq. The remaining jq = 1/2 states decay via S-waves and 
are expected to be quite broad. Although they have not yet been directly observed, their total production rate has 
been measured in i3-meson semileptonic decays |ic| . 

CLEO has observed the production of both of the narrow D** mesons in S ^ D*7nr decays with the following 
branching fractions [l^ : 

B{B- ^ D"^Tr-)xB{Dl^ D*+Tr-) = (7.8 ± 1.9) x 10"^, 
B{B^ ^ D^^TT-) X B{D;° ^ D*+7r-) = (4.2 ± 1.7) x 10"^ (1) 



The ratio of the B meson branching fractions 



is calculated in HQET and the factorization approach in Refs. [I3, |19|- In Ref. [18| R is found to depend on the 
values of sub-leading Isgur-Wise functions (fi_2) describing h-QCD/rric corrections. Variations of fi,2 by ±0.75 GeV 



B{B- DfiT-^ 
B{B- D^Ti-) 



result in values of R that ranee from to 1.5. In Ref. some of the sub-leading terms are estimated and the ratio 
is determined to be 



0.35 



1 + S, 



D2 

(3) 



Si" 



where S^^^^'^'' are non-factorized corrections that are expected to be small. The value of R calculated from the CLEO 



results given in Eq.lP) plus the ratio of branching fractions B{D^" ^ D+TT-)/B{Df D*+n-) ^ 2.3 ± O.i 
and the assumption that Di and D2 decays are saturated by the two-body Dn, D*Tr modes, is i? = 1.8 ± 0.8. This 
value is higher than the prediction, although the uncertainties are large. If more precise measurements do not indicate 
lower values of R, a problem for theory may arise. Thus, a measurement of R will allow us to test HQET predictions. 

Another possible inconsistency between theory and experiment is in the ratio of the production rates of narrow 
and broad states in semileptonic B decays. QCD sum rules predict the dominance of narrow D**{jq = 3/2) state 
production in B D**lv decays. On the other hand, the total branching fraction B{B D'^*^'n:l^v) — (2.6 ± 0.5)% 
measured by ALEPH and DELPHI Q is not saturated by the contribution of the narrow resonances, (G.86±0.37)% |2J , 
indicating a large contribution of broad or nonresonant D^*^7r structures. 

In this study we concentrate on charged B decays to £''*'^7r^7r^. For these decays the final state contains two 
pions of the same sign that do not form any bound states, making analysis of the final state simpler. 



THE BELLE DETECTOR 

The Belle detector [22 is a large-solid-angle magnetic spectrometer that consists of a three-layer silicon vertex 
detector (SVD), a 50-layer central drift chamber (CDC) for charged particle tracking and specific ionization measure- 
ment (dE/dx), an array of aerogel threshold Cerenkov counters (ACC), time-of- flight scintillation counters (TOP), 
and an array of 8736 CsI(Tl) crystals for electromagnetic calorimetry (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux return located outside the coil is instrumented to detect Kl 
mesons and identify muons (KLM). We use a GEANT-based Monte Carlo (MC) simulation to model the response of 
the detector and determine the acceptance 13 ■ 

Separation of kaons and pions is accomplished by combining the responses of the ACC and the TOP with dE/dx 
measurements in the CDC to form a likelihood C{h) where h — (vr) or (K). Charged particles are identified as pions 
or kaons using the likelihood ratio (PID): 

pm(K) = , , , ; PID(7r) ^ , ^^""^ , , = 1 - V\V>{K). 

At large momenta (>2.5 GeV/c) only the ACC and dE/dx are used since the TOP provides no significant separation 
of kaons and pions. Electron identification is based on a combination of dE/dx measurements, the ACC responses and 
the position, shape and total energy deposition {E/p) of the shower detected in the ECL. A more detailed description 
of the Belle particle identification can be found in ref. 



EVENT SELECTION 



A 60.4 fb~^ data sample (65.4 million BB events) collected at the T(4S') resonance with the Belle detector is used. 
Candidate B~ — > D+tt^tt^ and B^ _D*+7r^7r^ events as well as charge conjugate combinations are selected. The 
and D*^ mesons are reconstructed in the Z?+ K^tt^tt^ and D^tt^ modes, respectively. The D'^ 
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candidates are reconstructed in the D*^ K~tt~^ and D*^ K~tt^tt^tt~ channels. The signal-to-noise ratios for 
other D decay modes are found to be much lower and they are not used in this analysis. 

Charged tracks are selected with requirements based on the average hit residuals and impact parameters relative 
to the interaction point. We also require that the polar angle of each track be within the angular range of 17° — 150° 
and that the transverse track momentum be greater than 50 MeV/c for kaons and 25 MeV/c for pions. 

Charged kaon candidates are selected with the requirement PID{K) > 0.6. This has an efficiency of 90% for kaons 
and a pion misidentification probability of 10%. For pions the requirement PID(7r) > 0.2 is used. All tracks that are 
positively identified as electrons are rejected. 

D"*" mesons are reconstructed from K^n^n^ combinations with invariant mass within 13 MeV/c^ of the nominal 
mass, which corresponds to about SaK-mv- For mesons, the Kw or Ktttttt invariant mass is required to be 
within 15 MeV/c^ of the nominal mass {3aKTr)- We reconstruct D*+ mesons from the Dtt combinations with a 
mass difference of Af^,^ — M^o within 1.5 MeV/c^ of its nominal value. 

Candidate events are identified by their center of mass (cm.) energy difference, AE ~ ~ ^b, and beam- 

constrained mass, Afbc = — (X^jPO^i where E\, — is the beam energy in the T(4S') cm. frame, and pi 

and Ei are the cm. three-momenta and energies of the B meson candidate decay products. We select events with 
A/bc > 5.20 GeV/c^ and |A£;| < 0.10 GeV. 

To suppress the large continuum background {e^e~ qq, where q = u,d, s, c), topological variables are used. Since 
the produced B mesons are almost at rest in the cm. frame, the angles of the decay products of the two B mesons 
are uncorrelated and the tracks tend to be isotropic while continuum qq events tend to have a two-jet structure. We 
use the angle between the thrust axis of the B candidate and that of the rest of the event (Qthrust) to discriminate 
between these two cases. The distribution of | cos Qthrust\ is strongly peaked near | cos Qthrust\ = 1 for qq events and 
is nearly flat for T(45) BB events. We require | cosQthrust\ < 0.8, which eliminates about 83% of the contimmm 
background while retaining about 80% of signal events. 

There are events for which two or more combinations pass all the selection criteria. According to a MC simulation, 
this occurs primarily because of the misreconstruction of a low momentum pion from the D** D^*\ decay. To 
avoid multiple entries, the combination that has the minimum difference of Z coordinates at the interaction point, 
\Zt;-^ — ^^TTal, of the tracks corresponding to the pions from B D**t:i and D** D^*'^TT2 decays is selected 25]. This 
selection suppresses the combinations that include pions from Ks decays. In the case of multiple D combinations, 
the one with invariant mass closest to the nominal value is selected. 

B- D+Tv'TT- ANALYSIS 

The Afbc and AE distributions for B^ D+tt^tt^ events are shown in Fig. |21 The distributions are plotted for 
events that satisfy the selection criteria for the other variable: i.e. \AE\ < 25 MeV and |Afbc — Mb\ < 6 MeV/c^ for 
the A^bc and AE histograms, respectively. A clear signal is evident in both distributions. The signal yield is obtained 
by fitting the AE distribution to the sum of two Gaussians with the same mean for the signal and a linear function 
for background. The widths and the relative normalization of the two Gaussians are fixed at values obtained from 
the MC simulation while the signal normalization as well as the constant term and slope of the background linear 
function are treated as free parameters. 

The signal yield is 1101 ± 46 events. The detection efficiency of (18.2 ± 0.2)% is determined from a MC simulation 
that uses a Dalitz plot distribution that is generated according to the model described in the next section. 
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M^^ (GeV/c^) AE (GeV) 

FIG. 2: The a) Mhc and b) AE distributions for D^tx^-k~ events. The hatched histogram in (b) is the D mass sideband 

(||Mo -Mk^^I - 26MeV/c2| < e.SMeV/c^). 

Using the branching fraction B{D^ K~tt^h^) = (9.1 ± 0.6)% 0, we obtain 

B{B- D+TT-TT-) = (1.02 ± 0.04 ± 0.15) x 10~^ 

which is consistent with the upper limit obtained by CLEO, B{B^ D^tt^tt^) < 1.4 x 10"'^ The statistical 

significance of the signal is greater than 25a This is the first observation of this decay mode. The second 

error is systematic and is dominated by a 10% uncertainty in the track reconstruction (a 2% per track uncertainty 
was determined by comparing the signals for rj —>■ tt+tt^tt*' and r] —^ 77). The uncertainty in the _D+ —>■ A'^7r+7r+ 
branching fraction is 6.6% and that for the particle identification efficiency is 5%. Other contributions are smaller. 
The uncertainty in the background shape is estimated by adding higher order polynomial terms to the fitting function, 
which results in less than a 5% change in the branching fraction. The MC simulation uncertainty is estimated to 
be 3%. The possible contribution from charmless B-meson decay modes is estimated from the Mu sidebands. The 
sideband distribution, shown as the hatched histogram in Fig.E^b), indicates no excess from such events in the signal 
region. 



B — > Dm Dalitz plot analysis 

For a three-body decay of a spin zero particle, two variables are required to describe the decay kinematics; we use 
the two Dtt invariant masses. Since there are two identical pions in the final state, we separate the pairs with maximal 
and minimal Af^i^ values. 

To analyze the dynamics of i? ^ Dtttt decays, events with AE and Mbc within the |A£'| < 25 MeV, |Mbc — Mb\ < 
6 MeV/c^ signal region are selected. To model the contribution and shape of the background, we use a sideband 
region defined as 100 MeV > \AE\ > 30 McV with the signal Mbc given above. The minimal Dtt mass distributions 
for the signal and sideband events are shown in Fig. 13 where narrow and broad resonances are visible. 

The distributions of events in the A/|)^ versus M|,^ Dalitz plot for the signal and sideband regions are 
shown in Fig. ^ The Dalitz plot boundary is determined by the decay kinematics and the masses of their daughter 
particles. In order to have the same Dalitz plot boundary for events in both signal and sideband regions, mass- 
constrained fits of Kttt: to Md and ZJtttt to Mb are performed. The mass-constrained fits also reduces the smearing 
from detector resolution. 
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FIG. 3: a)The minimal D-k mass distribution of D^tv^tt^ candidates. The points with error bars correspond to the 

signal box events, while the hatched histogram shows the background obtained from the sidebands. The open histogram is 
the result of a fit while the dashed one shows the fit function in the case when the narrow resonance amplitude is set to zero. 
b)The background-subtracted Dtt mass distribution. The points with error bars correspond to the signal box events, hatched 
histograms show different contributions, the open histogram shows the coherent sum of all contributions. 



To extract the amplitudes and phases of different intermediate states, an unbinned fit to the Dalitz plot is performed 
using a method similar to CLEO's The event density function in the Dalitz plot is the sum of the signal and 
background. 

Since the Dtt mass distributions for the upper and lower halves of the AE sideband have similar shapes, we can 
expect similar background behavior for the signal and sideband regions. The background shape is obtained from an 
unbinned fit of the sideband distribution to a smooth two-dimensional function. The number of background events 
in the signal region is scaled according to the relative areas of the signal and the sideband regions. 

In the D+vr^vr^ final state a combination of the D^-meson and a pion can form either a tensor meson Z?2*' or 
a scalar state Dq^; the axial vector mesons D° and Df cannot decay to two pseudoscalars because of the angular 
momentum and parity conservation. The D*° cannot decay to the D+tt" because the -D*" mass is lower than that of 
D^TT^ . However, a virtual D* of a higher mass can decay to D+tt^. Another virtual hadron that can be produced in 
this combination is B* {B B*ti and B* Dtt). The contributions of these intermediate states are included in the 
signal-event density {S{q1, )) parameterization as a coherent sum of the corresponding amplitudes together with a 
possible constant amplitude (03): 

Siqlql) = |a^.A(2)(g^g|) +a^.e^*-SA(«)(g^g|) +a^,e^*-i +aB;^^^^^ 

(4) 

where ^TZ{Aq'^) denotes convolution with the experimental resolution. Each resonance is described by a relativistic 
Breit-Wigner with a q^ dependent width and an angular dependence that corresponds to the spins of the intermediate 
and final state particles: 

A'-'(4.4) = 4;g-(P.) ,;-MnX^,(,;, ^!o(P.) + (« - «). (5) 
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FIG. 4; The Dalitz plot for a) signal events and b) sideband events. 



where 

is the dependent width of the D** , with mass Ml and width F^, decaying to Dir state with orbital angular 
momentum L. The variables pi, p2, Pd, Qi = P2 +Pd, 12 = Pi + Pd arc the four-momenta of the pions, D, and 
Dtt combinations, respectively; P2, P2 are the magnitude of the pion three-momentum in the D** rest frame when 
the D** has a four- momentum-square equal to and M^, respectively; pi, p^ are the magnitude of the pion three- 
momentum in the B rest frame for the case when the D** four- momentum squared is equal to q^ and M£, respectively. 
The angular dependence for different spins of the intermediate states is: 

r(«)(gi,g,) = l, T(^)iq„q,) = M^,ose, r(^)(,i, = M|PiPi(cos^ - 1/3), (7) 

where 9 is the angle between the first pion from the B-decay and the pion from the D**-decay in the D** rest 
frame, and -F^'j^,, (pi) and J^^'Lj5(p2) are transition form factors, which are the most uncertain part of the resonance 
description. For the B D** and D** D form factors, we use the Blatt-Weiskopf parameterization 



where r=1.6 (GeV/c)^^ is a hadron scale. For the virtual mesons D* and B* that are produced beyond the peak 
region, another form factor parameterization has been used: 

i^AB(p) = e-'-(P-P°); (9) 

this provides stronger suppression of the numerator in Eq. ^ far from the resonance region. The resolution function 
is obtained from MC simulation; the detector resolution for the Dn invariant mass is about 4 MeV/c^. 

The D** resonance parameters {M^to, F^i.o, M^i.o, F^j.o) as well as the amplitudes for the intermediate states 
and relative phases (a£)» , ao', an*, as*, 03, 0_d*, 4>B', 03 ) are treated as free parameters in the fit. 
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TABLE L Fit results for different models. The model used to obtain the results includes amplitudes for D2, Dq,DI, 
intermediate resonances. Adding the constant term (ph.sp(a3)) does not significantly improve the likelihood. 

TablclII gives the results of the fit for different models. When the D* amphtude is included, the likelihood significantly 
improves and gives branching fractions values that are consistent with expectation based on the D* vi^idth and the 
D*'^TT^ branching fraction. When the B* amplitude is added, the likelihood is also significantly improved. A 
constant phase space term, 03 exp(i(^3), does not substantially change the likelihood and the final results are presented 
without this term. The variation of the fit parameters when these last amplitudes are included is used as an estimate 
of the model error. 

Figure 121 shows the Dtt mass distributions for different Dtt helicity regions. The helicity {cos 9h) is defined as the 
cosine of the angle between the pions from the B and D** decays in the rest frame of D** . The number of events in 
each bin is corrected for the MC-determined efficiency. The curve shows the fit function for the case when the D2 , Dq , 
D* and B* amplitudes are included. The D2 resonance is clearly seen in the helicity range |cosf?/i| > 0.67, where 
the D-wave peaks. The range 0.33 < | cos6'/i| < 0.67 where the D-wave amplitude is suppressed, shows the S-wave 
contribution from the Dq while the low helicity range | cos 0ft, | < 0.33 demonstrates a clear interference pattern. 

Another demonstration of the agreement between the fitting function and the data is given in Fig. |H| where the 
helicity distributions for different regions are shown. The histogram in the region of the D2 meson clearly indicates 
a I cos^ 6'/,. — 1/3 p D-wave dependence. The distributions in the other regions show reasonable agreement between the 
fitting function and the data except for a few bins in the small Mu-^ min region and with helicity close to 1 (Fig.j^a)). 
This region is populated mainly by the virtual D* and B* production, the description of which depends on the form 
factor behavior. This discrepancy does not affect the determination of the D** parameters that are the main topic of 
this work. 

The fit quality is estimated using a two-dimensional histogram of minimum versus the Di: helicity and calcu- 
lating the /N for the function obtained from unbinned likelihood minimization. The confidence level for the model 
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FIG. 5: The minimal Dtt mass distribution for different helicity ranges. The two curves are the fit results for the case of 
D2, Dq and D* amplitudes (the top curve) and the background contribution (the bottom one). The number of events in each 
bin is corrected for the efficiency obtained from MC simulation. 



Model 


-2ln{C/£max) 


D2 I Dq , D-u , By 





D2, Dl, B*, ph.sp(a3) 


265 


D2, Dl , B*, 1" 


355 


D*2, Dl , B*,2+ 


235 



TABLE II: Comparison of models with and without a O"*" resonance. The amplitudes for D| and the virtual Dy and By are 
always included. 



with D'2, Dq, Z?* and B* is about 0.8%. The low confidence level is due to the poor description in the region where 
Mdtt min is Small and Mo-n- max is large (or helicity is close to 1) as discussed above. 

In Table Hll the likelihood values are presented for the case when the broad scalar resonance is excluded or when it 
has quantum numbers different from — 0+. For all cases the likelihood values are significantly worse. Thus, we 
claim the observation of a broad state that can be interpreted as the scalar Dq . The fit gives the following parameter 
values: 

M^i-o = (2308 ± 17 ± 15 ± 28)MeV/c^ T^.o = (276 ± 21 ± 18 ± 60)MeV. 

The values corresponds to the case when four amplitudes (column III of Tabled) were included. Here and throughout 
the paper the first error is statistical, the second is systematic and the third is the model-dependent error described 
below. 

The values of the narrow resonance mass and width obtained from the fit are: 

Af^.o = (2461.6 ± 2.1 ± 0.5 ± 3.3)MeV/c2, T^'O = (45.6 ± 4.4 ± 6.5 ± 1.6)MeV. 
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The value of the D'^ width is larger than the world average of 23 ± 5 MeV, and is consistent with the preliminary 
result from FOCUS of 30.5 ± 4.2 MeV The previous analyses did not take the interference of intermediate states 
into account; this suggests that there may be large unaccounted systematic errors in these measurements. 
The following branching ratio products are obtained: 

BiB-^ -> Dfn-) X B{Df D+TT-) = (3.4 ± 0.3 ± 0.6 ± 0.4) x 10"'', 

B{B- D*^TT-) X B{d;" D+TT-) = (6.1 ± 0.6 ± 0.9 ± 1.6) x 10"^, 

and the relative phase of the scalar and tensor amplitude is 

0^.0 = -2.37 ± 0.11 ± 0.08 ± 0.10. 

The systematic errors are estimated by comparing the fit results for the case when the background shape is taken 
separately from the lower or upper sideband in the AE distribution. The fit is also performed with more restrictive 
cuts on AE, Mhc and AMd. The maximum difference is taken as an additional estimate of the systematic uncertainty. 
For branching fractions, the systematic errors also include uncertainties in track reconstruction and FID efficiency, as 
well as the error in the K~-k'^'k^ absolute branching fraction. 

The model uncertainties are estimated by comparing fit results for the case of different models (II-IV in Tabled 
and for values of r that range from to 5 (GeV/c)^^ for the transition form factor defined in Eqs. ((H)) and 

B- D'+Tf-TT" ANALYSIS 

For D* reconstruction, the _D*+ D°7r+ decay is used and two decay modes K~tt^ and iir~7r+7r+7r~ 

are included. The AE and Afbc distributions are shown in Fig. [3 In each mode the number of signal events is 
obtained in a way similar to that described for the Dtttt selection. The observed signal yields of Nktt = 273 ± 21 and 
Nk3tt = 287 ± 22 for the Ktt and Kmnr modes, respectively, are consistent, based on the D branching fractions and 
the efficiencies determined from MC: (13.6 ± 0.2)% for K-t:+ and (6.5 ± 0.2)% for K-tt+tt+t:-. 

The branching fraction of {D* — + 011)7111 events, calculated from the weighted average of the values obtained for 
the two modes, is: 

B{B- D*+Ti-TT-) = (1.25 ± 0.08 ± 0.22) x lO'^, 

where the first error is statistical and the second is systematic. This measurement is consistent with the world average 
value (2.1±0.6) x 10~^ 1]. The systematic error is dominated by the uncertainties in the track reconstruction efficiency 
(16 %) (for a low momentum track from the D* decay the efficiency uncertainty is 8 %) and the FID efficiency (5 %). 
The background shape uncertainty is estimated in the same way as for the ZJtttt analysis to be 5%. 

B D*iTn coherent amplitude analysis 

In this final state we have a decaying vector D* particle. Assuming the width of the D* to be negligible, there are 
two additional degrees of freedom and, in addition to two D*i: invariant masses, two other variables are needed to 
specify the final state. The variables are chosen to be the angle a between the pions from the D** and D* decay in 
the D* rest frame, and the azimuthal angle 7 of the pion from the D* relative to the B D*nn decay plane. 
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FIG. 6: The helicity distribution for data (points with error bars) and for MC simulation (open histogram). The hatched 
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below D2 resonance, c) - region of the tensor resonance, d) - region higher of the D2. 
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For further analysis, events satisfying the selection criteria described in the first section and having AE and Afbc 
within the |A£;| < 30 MeV, |Mbc - Mb\ < 6 MeV/c^ signal range are selected. To understand the contribution and 
shape of the background, we use events in the 100 MeV > \AE\ > 30 MeV sideband. 

The D*Tr final state can include contributions from the narrow £'2° and Dj, and the broad D'^ states. The 



13 



minimal D*7r mass distributions for the signal and sideband events are shown in Fig.|Sl A narrow structure around 
Mjj'tt ^ 2.4 GeV/c^ and a broader component that can be interpreted as the D[ are evident. 

The Dalitz plot distributions for the signal and sideband events are shown in Fig. |^ In order to have the same 
boundary of the Dalitz plot distributions for events from both signal and sideband regions as well as to decrease the 
smearing effect introduced by the detector resolution, mass-constrained fits of Dn to Mjj,+ and D*TnT to Mb are 
performed. 




2.2 2.3 2.4 2.5 2.6 2.7 2.2 2.3 2.4 2.5 2.6 2.7 

FIG. 8: a) The minimal mass distribution of D*n events. The points with error bars are experimental data, the hatched 
histogram is the background distribution obtained from the sideband, the open histogram is MC simulation with the amplitudes 
and intermediate resonance parameters obtained from the fit. The dashed histogram shows the contribution of the broad 
resonance, b) The background-subtracted D*Tr mass distribution. The points with error bars correspond to the signal box 
events, hatched histograms show different contributions, the open histogram is a coherent sum of all contributions. 



To extract the amplitudes and phases for different intermediate states, an unbinncd likelihood fit in the four- 
dimensional phase space was performed. Assuming that the background distribution (^(g^, gf, a, 7)) in the signal 
region has the same shape as in the AE sideband, we obtain the B{ql,q2,a,j) dependence from a fit of the side- 
band distribution to a smooth four-dimensional function. The number of background events in the signal region is 
normalized according to the relative areas of the signal and the sideband regions. The signal is parameterized as a 
sum of the amplitudes of an intermediate tensor {D2), and two axial vector mesons (-D'l, -Di) convoluted with the 
resolution function Tl{Aq^) obtained from MC simulation: 

S{ql,qla,j) = \A^^^\ql qj, a,^) + A^'^'^^ql qj, a,j) + A^'^'^Hql ql a,j) (10) 
+ ai3„ e'^°^- A^^' {ql , <Z2^ «, 7) + as; e'^«.* A^* {qf , <Z2^ «, 7) + a^e'^' \ ' ® R{Aq^). 

Each resonance is described by a relativistic Breit-Wigner with a width depending on q^ (see Eq|S)l and angular 
dependence corresponding to the spins of the intermediate and final state particles. 

Al"^ p^prP3 

T''^^\qi,q2,ce,j) = ajj^ — ^ ^ (sin cos 7 sin a -|- 2 cos cos a), 
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FIG. 9; The Dalitz plot for a) signal and b) sideband events. 



1^ '[qi,q2,a,^) = a^ee i ■= — (sm 6^ cos 7 sm a — cos y cos a j 



(11) 



r(^^*) (91,92, a, 7) 



i4>o'0 MbpIpIp3 



cos a sm tl sm a sm 7, 



where a^ij, a^i^, a^)*, 4'd[^ (t>D* sue the relative amplitudes and phases for transitions via the corresponding inter- 
mediate state. The amplitudes of S and D waves in Eq. pi|l correspond to decay via 1^^^ ^-nd ly^ intermediate 
states, respectively. Due to the finite c-quark mass, the observed 1+ states can be a mixture of pure states. Thus, 
the resulting amplitude will include a superposition of the amplitudes for the corresponding Breit-Wigner: 



r(^i)(9i, <72, a, 7) - r(i^)(gi, 92, a, 7) cos w - e'^'T^^''\quq2, a, 7) sinw 
r(^i)(9i, 92, a, 7) = r(i^)(9i, 92, a, 7) smco + e~'^T^'''\qi, q^, a, 7) cos 



(12) 



where uj is the mixing angle and ip is a complex phase. 

The D*TT pair in the final state can be produced via virtual or B*" decaying to D*^tt^ . Inclusion of a virtual 
significantly improves the likelihood; including in addition B* and a constant term also improved the likelihood, but 
the significance is not high (see Table ^Oj. A fit without the inclusion of a broad resonance gives a considerably worse 
likelihood (see Table irv|) . We also tried to fit the data by including a broad resonance with other quantum numbers 
such as 0~, 1~, 2+. In these cases the likelihood is also significantly worse, as shown in Table HVl We conclude that 
we have observed the broad 1+ D[ state with a statistical significance of more than 10 a. The model and systematic 
errors are estimated in the same way as for the Dnn case. 

The D2 mass and width are fixed to the values obtained from the Dmr analysis. The axial vector D** masses 
and widths as well as the branching fractions and phases of amplitudes a/j^, a^jj, cld*, (t'D[i 4'd* are treated as free 
parameters of the fit as are the mixing angle ld and the mixing phase "0. 

Since there is no good way to graphically present the data and the model in four dimensions, we show the projections 
of the distributions for various variables. Figure|Hlshows the Md'-k mm distribution together with MC results that were 
generated according to the model containing Di , D'l , D2 and virtual , B* intermediate resonances with parameters 
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TABLE III: Fit results for different models. The model that is used to obtain these results includes amplitudes for 
D2 , Di , D'l , Dv , By intermediate resonances. Adding a constant term does not improve the likelihood significantly. 

obtained from the fit. Figure [TUl shows a comparison of the data and the MC simulation for D** and D* hehcities as 
well as the angle 7 for ranges corresponding to the two narrow resonances Di {q^ — (5.76 ~ 5.98) (GeV/c^)^), D2 
[q^ = (5.98 ^ 6.15) (GeV/c^)^) and the regions populated mainly by the broad D[ state below [q^ < 5.76 (GeV/c^)^) 
and above {q^ > 6.15 (GeV/c^)^) the narrow resonances. All distributions indicate good agreement between the data 
and the fit result. We cannot characterize the quahty of the fit by the standard test since for a binned distribution 
with four degrees of freedom and a limited data sample any reasonable binning will result in only a few events per 
bin. Therefore to estimate the quality of the fit we determine values for different projections of the distributions 
in Figs. island 1101 The obtained values correspond to confidence levels in the 5 — 90% range. 
For the Di meson we obtain the following parameters: 

M^o = {2421 A ± 1.5 ± 0.4 ± 0.8) MeV/c^ T^,!. = (23.7 ± 2.7 ± 0.2 ± 4.0) MeV. 

These parameters are in good agreement with the world average values: AIjjo — (2422.2 ± 1.8) MeV/c^, Tjjo — 
(18.9j:^J)MevQ. 

The broad D'^ resonance parameters are: 

Mjj.o = (2427 ± 26 ± 20 ± 15) MeV/c^ F^,™ = (384j:^f ± 24 ± 70) MeV. 

Observation of a similar state was reported by CLEG but was not published; our measurement is consistent with 
CLEG'S preliminary results: Mjjm = (2461^^^) MeV/c^, T^'o = (290t^J°) MeV 0- 
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TABLE IV: Comparison of the models with and without a broad 1"*" resonance. The D2 and Di amphtudes are always included. 

The results for the products of the branching fractions of the B and D** mesons are: 

B{B- Dnr-) X B{Di D*+7t-) = (6.8 ± 0.7 ± 1.3 ± 0.3) x 10"*, 

B{B- Dfn-) X B{Df D*+tt-) = (1.8 ± 0.3 ± 0.3 ± 0.2) x 10""*, 

S(B^ ^ D[°7r-) X 6(L'f ^ D*+7r-) = (5.0 ± 0.4 ± 1.0 ± 0.4) x lO"'', 
the relative phases of the and D[ amplitudes are: 

0Df = -0-57 ± 0.14 ± 0.06 ± 0.13; (f>D[ = 1.68 ± 0.20 ± 0.07 ± 0.16, 
and the mixing angle of two axial states and the complex phase are: 

cj = -0.10 ± 0.03 ± 0.02 ± 0.02, 

= 0.05 ± 0.20 ± 0.04 ± 0.06. 

To understand the uncertainties in the background shape and the efficiency of the cuts, additional studies are 
performed. The background shapes obtained separately for the upper and lower AE sidebands are used in the 
likelihood optimization. We also apply more restrictive cuts on AE and Afbc that improve the signal-to-background 
ratio by about a factor of two and repeat the fit. The maximum difference between the values obtained with different 
cuts and different background shapes is included in the systematic error. The branching fraction errors also include 
an 18% systematic uncertainty in the detection efficiency. The model uncertainties are estimated by comparing fit 
results for the case of different models (Table ^O} and for values of r in the range from to 5 (GeV/c)^^, where r is 
the hadron scale parameter in the transition form factors of Eqs. © and ©. 



DISCUSSION 

From the measured products of branching fractions of B{B^ Dj"^ )'S(-C2° ^ _D*+7r^) and B{B^ — > 
D2^Ti^)B{D2^ D+TT^) we obtain the ratio of the D^^ branching fractions: 

B{Df D*+TT-) 

which is consistent with the world average H = 2.3± 0.6. Theoretical models Q H Q 

predict H to be in the 

range from 1.5 to 3. If the D2 decay is saturated by Dtt, D*it transitions, and Di decay by D*ti, then the ratio R 
in Eq. ^ can be expressed as the following combination of branching fractions: 

^ = B{B-^D-,.-)B{D-,^D*+n-) ^ ^ 
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FIG. 10: Distributions of the helicity angle of D** (cosd), helicity angle of D*(cosQf) and azimuthal angle 7 in four different 
regions. The points are the experimental data, the histogram is MC simulation with fitted parameters, the hatched histogram 
shows the background contribution (from the AE sideband). 
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The obtained value is lower than that of the CLEO measurement (although the measurements are consistent within 
errors) but is still a factor of two larger than the factorization result [l^ . From our measurement it is impossible to 
determine whether the non-factorized part for tensor and axial mesons is large, or whether higher order corrections to 
the leading factorized terms should be taken into account. According to Ref. 18], the observed value of i? corresponds 
to a value of the sub-leading Isgur-Wise function fi ~ 0.401q ijg GeV. 

For semileptonic decays, where there is no non-factorized contribution, the corresponding ratio is 0.5 ± 0.6 |l2| . 
which, within experimental errors, is consistent with both our measurement and the model prediction. More accurate 
measurements of semileptonic modes containing D** mesons may help resolve this problem. 

Our measurements show that the narrow resonances comprise (36 ± 6) % of Dtttt decays and (63 ± 6) % of D*TnT 
decays. This result is inconsistent with the QCD sum rule [20j that predicts the dominance of the narrow states in 
B D'-*^Tnr decays. It is also possible that in i?" D**^tt~ decays the color suppressed amplitude is comparable 
to the tree amplitude, so that other transition form factors play a role. The ratio of the production rates for narrow 
and broad Dt: states in semileptonic B D^*^'kIi> decays measured at LEP 12] also indicates an excess of the broad 
states. More accurate measurements of both semileptonic decays and other charged states of the 0^*^7177 system may 
resolve this discrepancy. 



CONCLUSION 

We have performed a study of charged B D^tt^tt^ and B D*^TT^n^ decays. The total branching fractions 
have been measured to be B{B- D+tt^tt") = (1.02 ± 0.04 ± 0.15) x 10"^ and B{B- D*+TT-Tr-) = (1.25 ± 
0.08 ± 0.22) X lO^'^. For the former decay this is the first measurement. 

A study of the dynamics of these three-body decays is reported. The D^tt^tt^ final state is well described by the 
production of -DJtt^ and DqTt^ followed by D** Dir. From a Dalitz plot analysis we obtain the mass, width and 
product of the branching fractions for the -Dj": 

M^.o = (2461.6 ± 2.1 ± 0.5 ± 3.3) MeV/c^ F^.o = (45.6 ± 4.4 ± 6.5 ± 1.6) MeV, 

B{B- -> D;°Tr-) X B{D;" D+TT-) = (3.4 ± 0.3 ± 0.6 ± 0.4) x 10"''. 
In this mode we also observe production of a broad scalar Dq meson with mass and width: 

M^.o = (2308 ± 17 ± 15 ± 28) MeV/c^ Tjj.o = (276 ± 21 ± 18 ± 60) MeV. 
The product of the branching fractions for the Dq state is 

S(B' ^ D^°7:-) X B{Dl° ^ D+TI-) = (6.1 ± 0.6 ± 0.9 ± 1.6) x 10"'', 
and the relative phase of the scalar and tensor amplitudes is: 

,/,^.o = -2.37 ± 0.11 ± 0.08 ± 0.10. 

This is the first observation of the D^. 

The D*TnT final state is described by the production of D2Tt, D[tt and Ditt with D** D*Tr. From a coherent 
amplitude analysis we obtain the mass, width and product of the branching fractions for the Di: 

Mdo = (2421.4 ± 1.5 ± 0.4 ± 0.8) MeV/c^ V = (23.7 ± 2.7 ± 0.2 ± 4.0) MeV, 
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B{B- DiTT-) X B{Di D*+-K-) = (6.8 ± 0.7 ± 1.3 ± 0.3) x 10"^, 
and measure the product of the branching fractions for the tensor meson process: 

B{B- D^°TT-) X B{D*2° D*+Tr-) = (1.8 ± 0.3 ± 0.3 ± 0.2) x 10-^ 
and the relative phase of the tensor meson to the axial vector Dj: 

(j)jj,o = -0.57 ± 0.14 ± 0.06 ± 0.13. 
We observe the broad D[ resonance with mass and width: 

Mjym = (2427 ± 26 ± 20 ± 15) MeV/c^ Tjj,o = (384l^f ± 24 ± 70) MeV. 
The product of the branching fractions is: 

B{B- D[\-) X B{D'° D*+n-) = (5.0 ± 0.4 ± 1.0 ± 0.4) x 10"^ 
and the relative phase of D'^ to D\ is: 

(j)D'^ = 1.68 ± 0.20 ± 0.07 ± 0.16. 
Our analysis also indicates that the axial vector states mix. The mixing angle is 

uj = -0.10 ± 0.03 ± 0.02 ± 0.02, 

and the phase is 

ij; = 0.05 ± 0.20 ± 0.04 ± 0.06. 
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